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ABSTRACT 


A pneumatic, pure fluid control system, capable of main- 
taining the angular position of an inertia load in a gravita- 
tional space, using Only one moving part, was designed and 
built. The control system consisted of a vortex-sink angular 
rate sensor, a pendulum controlled flapper-nozzle angular posi- 
tion sensor, a summing amplifier and bistable reaction jets. 
The system was designed to operate in a bang-bang, manner in 
order to reduce errors in a minimum amount of time. 


A previous author's work on the vortex-sink angular rate 
sensor is extended to include information and date on propor- 
tional amplification of the output signal and the frequency 
response of the sensor. Methods for synthesizing the components 
into a working system are presented along with a theoretical 
analysis of the system's performance. 


Results of experiments performed on the completed system 
showed that the reaction jets could be switched in a nearly op- 
timal manner when responding to errors introduced on one side 
of the system, but, due to the unsymmetrical gain of the rate 
sensor, the errors introduced on the opposite side resulted in 
a certain amount of after end-point chattering. 


Good correlation between theory and experiment led to the 
conclusion that the application of standard control system design 
techniques can be readily applied to the design of pure fluid con- 
trol systems. 


Thesis Supervisor: Herbert H. Richardson 
Title: Associate Professor of Mechanical Engineering 
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CEAPTER 1 


INTRODUCTION 


1.1 General Introduction to Pure Fluid state Technology 


The basico concept of pure fluid amplification is that 
high energy fluid streams can be controlled by low energy 
Streams. The concept was evolved at the Harry Diamond Lab- 
oratories of the U. S. Army Material Command, Washington, 
ECC. in 1959. $ 

At the First Fluid Amplification Symposium held at 
the Laboratories in October of 1962, Brigadier General 
John G. Zierdt described the possible future of these compon- 


ents: 


This aavance has removed many of 
the limitetions on conventional 
military and commercial pneumatic 
systems and has enabled an extrene- 
ly rapid extension of their capabil- 
ities. The use of low energy flow 
to direct and control high energy 
Streams permits the design of pure 
pneumatic components, elements, and 
Systems, capable of amplification, 
feedback, memory, logic functions, 
analogue computation, anc digital 
computation. This, in turn, makes 
it possible to use. pneumatic control 
systems, amplifiers, and computers 
where they could not be used before. 


1.2 Purpose anda Scope of Thesis : 


a 


Since the inception of pure fluid amplification, much 





work has been done in analyzing individual com»onents.  ite- 
cently the emphasis in the field has been turned to the ĉe- 

| velopment of complete systems. Although much worx has been 
done along these lines, the quantitative data associated vith 
‚such systems is scanty and iena 

One objective of this thesis is to develop techniques 

that may be. used in designing, enalyzing and synthesizing, these 

components into a working control system. 

The specific objective is to develop a pneumatic con- 
trol system, using as few moving parts as possible, capable 
of maintsining a commanded angular position of a pure inertia 
load. Should errors be introduced, the control system is to 
reduce these errors to near zero in the shortest possible 
time within the limitations of the prime move used. 

Another design objective is to purposely make the 
system respond slowly so that the device may be used as an 
educational model demonstrating the working of individual 


components and the system as a whole. 


1.5 Discussion of the Proposed System to attain the Design 
Cbjectives | 


Minimum time control of a vehicle generally involves 
the use of the full force available from the prime mover in 
a béng-bang manner. for example, it is obvious that to move 
an automobile from a point, A, to another point, B, in a min- 
imum amount of time, the operator should use maximum acceler- 


etion, at the right instant lock his brakes and slide into the 


A 4 





desired end point. The operator in this example would De 
bang Pull driving force in one direction ane full braking 
force in the other to achieve the minimum time ob jective. 

in controlling an anguler inertie load, it can be 
shown that minicun time movements sre obtained by using full 
torcue in one direction or the other. igure 1.1 shows a 
phase plane portrait of an inertia load being driveh by a 
constant torque, whose sense is either plus or minus, for 
various initial conditions. 

It can be seen from Figure 1.1 that one curve exists, 
which, if usec as a boundary to switch the sense of the torque, 
will allow the system to reduce any errors to zero in mini- 
mum time. In order to switch the sense of the applied torque 
at the proper instant, it is obvious that the control system 
must be able to 1) measure angular positions, 2) measure 
angular rete, 3) compute the coordinates of the switching 
curve, and 4) conpare the values of angular rate and anguler 
position with the values of the switching curve and switch 
the sense of the applied torque at the proper instant. 

| The remaining chapters of this thesis deal with the 
design of the components needed to accomplish the above men- 
tionec functions and the synthesis of these components into 


a working system. 


T Superscripts csuute references in the Bibliography 
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CHAPTER II 


VORTEX ANGUL.R RATE SENSOR ANALYSIS AND DESIGN 


2.1 General Theory of Operation 


The principles on whick a vortex engular rate sensor 
operate. are quite simple. Real difficulties quickly Sea 
however, when attempts are MES at mathematical analysis of 
the complex flow which takes place in these devices or when 
attempts are made to obtain pressure or flow signals which are 
some function of the input angular rotation. It is Dot 
because of the latter difficulty that the methods of obtain- 
ing output signals have been, and still are, held confidential. 

The basic operation of a rate sensor is depicted in 
Figures 2.laend 2.l1b. Figure 2.la shows the flow state ín the 
absence of angular rotation. In this case a particle of 
fluid passes from tbe supply plenum through the porous ring 
into the vortex chamber and out the exit hole. From mass 


continuity the radial flow velocity is described by 


Ur = mo S I (2.1) 
ar n 
where 


q = flow rate supplied to the sensor 


ut 
I 


- any arbitrary radius 


Do 
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thickness of the vortex chamber 
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Thus just inside the porous ring, the radial velocity 


is given by 


Uo = S fone) 


and at the exit 


Ue - Q e (233) 
. a h 


where Q and h are as given above 


radius of porous ring 


and Re 
Ro 


exit hole redius 

Figure 2.1b shows the flow condition when the sensor 
is rotated at a constant anguler rate. In this cese a par- 
ticle of fluid entering the vortex chanber has the radial 
velocity given by equation (2.1) and a tangential component 
of velocity imparted to it while passing through the porous 
ring. In visualizing this tangential component, one must 
use a frame of reference other than the vortex chamber. Thus 
with respect to laboratory coordiínetes, the magnitude of the 
tangential component of velocity is given by 

Uto = RoW (274) 
where Uto = magnitude of the tangential velocity 
W = magnituce of the input rotation 

Hence, a particle of fluid of mass dm will have angular mo- 
mentum given by Ray Uto dm. Conservation of Anguler Momentum 
produces an expression for the tangential velocity (Ut) at 


an: arbitrary radius (r): 


Tt Ut dm = Ro Uto dm (235a) 
Ur = _Ro2W (2.50) 
r 





For an ideal flow the farilisr result of U. aoproaching infin- 
ity is seen as r approaches zero. 

Thus it can be seen that if a means can be devised to meas- 
ure the amount of tangential velocity directly or indirectly at 
So axrticular radius, this' velocity will change in proportion te 
the magnitude of @as to tne direction of rotation. 

The method used to measure the amount of tangential veloc- 
ity present in the vortex chariber of the rate sensor used in this 
control system was developed by Sarpkaya*. The analytical analy- 
sis about to be discussed is also given by Sarpxaya but is included 
here to maintain continuity of this chapter and to set the back- 
ground for the extension of Sarpkaya's work which follows in the 
remaining articles”. 

A schematic of the sensor developed by Sarpkaya is shown 
in Figure 2.2. This sensor differs from the one shown in Figure 
2.1 in that two sink-tubes are attached to the exit holes of the 
vortex chamber. These sink-tubes serve the following purposes: 

1) Tangential flow amplification - each tube can now be 
considered as a miniature elongated vortex chamber being fed from 
the top by a uniform radial flow on which is superimposed the 
tengential swirl flow from the main vortex chamber, 

2) Directing radial flow - the tubes receive the flow 
from the main vortex chamber and direct it toward the probe, 

3) The sink-tubes provide a housing for the probe. 

The probe, shown in Figure 2.3, consists of a small cir- 
cular cylinder divided into two chambers. Two small holes are 


drilled perpendicular to the center line of the probe and 


* To further credit the work of Sarpkaya, most of his 


Symbols are maintained throughout this chapter. 
8 
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Figure 223 


Cross Section of Probe 
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1 


aligned along a straight line drawn on the surface of the 
cylinder between its ends. The probe is inserted into the 
sink-tube perpendicular to the axis of flow, es shown in 
Figure 2.2, The holes are set at an angle to the oncoming 
flow of fluid in the sink-tube. 

The operation of this probe in sensing angular eee 
tion is as follows: From potential flow theory the pressure 
distribution around a circular cylinder placed in a uniform 


flow field is given by 


2 
Pp = Po ¢ Pus (1 4sin*B)' (2.6) 
where Py 7 pressure in lbf/in^ 


Po = static pressure in Lof/in® 


density in lbf-sec? 


D 
' 


in4 
" Ug * axial sink velocity in in/sec 
g = angular position measured erouná the 


cylincer from the stagnation point. 
& sketch of equation (2.6) is shown in Figure 2.4. in the 
absence of angulex rotation of the device, the uniform radial 
velocity in tne vortex chamber (U) turns into the sink-tube 
(now redesignated Us), strikes the probes and then passes out 
of the sensor. DBecsuse both holes are aligned at the same 
angular position along the probe, equel pressures exist in 
each chamber., The output signal of the probe, however, is 


the difference between these two pressures. Since both 
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chambers are at the same pressura, no output signal (Ar) 
exists. 

Now let there exist angular rotation of the rate sensor. 
Figure 2.5 shows the velocity vectors in way of the probe. 
The tangential velocity (vpo) now present in the sink-tube 
has the effect at diverting the main axial sink flow. The 
top hole is thus placed in a uniform flow coming from an 
angle which is greater than that which existed without angu- 
ler rotation. the bottom hole iS placed in a flow coming from 
an angle which is Smaller. Thus, the pressure in the chamber, 
to which the top hole is connected, decreases, (point A in 
Figure 2.4), while the pressure in the bottom chamber increases, 
(point B in Figure 2.4), A pressure differential now exists 
across the probe which is a function of the magnitude of the 
angular rotation. It is readily visualized that had the angu- 
lar rotation been in the opposite direction that the reverse 
procedure would have taken place. Hence, the signal also 
senses direction of anguler rotation as well as magnitude. 
It is assumed in this explanation that the relative magnitude 
of the sink-tube axial and tangential velocities are such 
that a long spiraling flow exists at the radius rpo. 

lo visualize the flow quantitatively and to obtain the 
function relating Ap anc Wrefer again to Figure 2.5. By 


maopection it can be seen that AG, Us, and vpo are related 


by: 


^0- tan"! ever (2.7a) 
Ug 
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sink-tube 


+ 





if AQ is a small angle then: 


V 

nM (2.76) 
Us 

To obtain the raximum signal (f\.P) possible for a 

given W,it is necessary to cetermine the maximum slope of 


the curve shown in Figure 2.4. Hence, cifferentiating 


equation (2.6) to obtain the slope gives: 


om |» pus? 8 sinQCos Y (2.8) 
30 2 


Now the derivative of equation (2.8) set to zero will give 





the Ó for maximum slope. 


s = | 
Ə P, = .g PUs Cos 29, 9=T/4 
—— ? 
30 * 2 2 
| (2.9) 
Therefore, the probe holes should be located at a 
45 degree angle to the oncoming flow. This value substituted 


into equation (2.8) gives 


IP, = -4/ Pu,* (2.10) 
2o 


39 


and in differential form for one hore: 


Ae! = + (Pu ^8 (2.11) 


Equation (2.5b) gives the relation between the tangential 
velocity existing in the sink-tube anà the input angular 


rotation. ‘With new symbols, equation (2.5b) becomes: 


Yon 7 Ro ¿O Mean 





Combining equations (2.7b), ana(2.11), and (2.12) gives: 


laz] = 4 RUY  , for one hole (2.13a) 
U 
a Usl yo 

2 


Since two holes are being used equation (2.16a) becomes: 
LA p| = 8 R y (2.180) 
2 Eom 
Pu p 


2 


ST po 

This is the desired relationship between Á P and Win non- 
dimensional form. Note that for a given sensor and probe cone 
figuration, that Ro and Doo 


are fixed. Further, if ihe flow 
4 , 
through the sensor is fixed, then Ug is constant. Thus it is 


seen that equation (2.13b) is of the form: 


AP = KW (2.14) 
where X > 8 [omues Becr , 8 constant 
—  — O PON ARO 
E Usro geometry and 
flow. 


Data given by Sarpkaya bear out the linearity of this ecuation 
over a substantial range of QJ. AS might well be expectec, 
from inspection of Figure 2.4, linearity is preserved proviaed 
A O is small. 

Thus far the equations derived to describe the opera- 
tion of this rate sensor, are based on an idealized potential 
flow. lt is now appropriate to consider the effects of vis- 
cosity on the operation of the sensor and the describing equa- 
tions. Two ma jor effects of viscosity will be noticed: 


1) Viscosity will cause the vortex motion in the 


16 





vortex chamber and sink-tube to deviate from the ideal vortex 
given sy equation (2.5b). 

2) Viscosity will cause the flow around the probe to 
Separate resulting in an aiteration of the pressure distri- 
bution around US probe given by equstion Bee) | 
Effect 1) is studied in great detail by Sarpkaya, by analyzing 
the boundary layer equations in the chamber for the radial flow 
with smell tangential flows superimposed on the main flow. 

The results show that viscosity reauces the strength of circula- 


iones the flow travels from R. to [poe Sarpkaya introduces 


a factor he cells the viscous efficiency of the sensor defined 


as; 
p =n (2.15) 
Po 
where È = strength of vortex circulation at 
radius r 
[s = strength of vortex circulation at 


racius Ry 


The viscous efficiency between the porous ring and the probe 


is: 
lo Rol 2rr Ro Ro° 
(2.168) 
from which 
= l 
DM QE NE (2.16b) 
po 
Poo 
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With equation (2.16b) used in lieu of equation (2.12) and then 
combined with esurtions (2.7b) and (2.11), the new relation- 


ship between AP and (y becomes 





pa z 8 ag 4022 for two holés. (2.17) 
P 3 Us Tro 
2 


Intuitively one might expeot this efficiency to be relatively 
high provided the circulation is quite small. This is shown to 
be true by Sarpkaya in a plot of E versus Reynolds a 
Effect 2) is an extremely important consideration which 
seems to have been neglected by Sarpkaya in his analysis of 
the sensor. Non idesl flow around the probe radically re- 
duces the slope of AP versus O in the vicinity of g= 45 
For the case of ideal flow arounà the probe, the nondi- 
mensional slope is given by equation (2.11). 
[Ar] / Lu _ o 
5 7 ES co 
A 6 | repeated 
From actual measurements taken on the probe used in 
the sensor constructed by the author, this slopb wes found to 
vary between 2.58 and 3.17 for the flows considered. The ef- 
fect of viscosity is more clearly shown in Figure 2.6 where 
nondimensional pressure distribution for two actual flows 


versus D is shown along with that predicted by potential 


flow theory. Thus the relationship between AP and W with 


the introduction of the new slope of Ap vs Ab Becomes: 


+ 
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5 : 
I^» | = 9.16 BR for two holes. 
Rug*/2 Us Typo (2,18) 


2.2 Rate Sensor Design and Desizn Procedure* 


The author did not consider the effect of slope reduce 
tion due to viscosity of the AP versus Q curve, Giscussed 
at the end of article 2.1, until after the rate sensor haá 


been designed. Hence, the equations used for the design were 


2 
f Z p de Pus AG: two holes (2.19) 
2 
| Ap = Ro? (QE (2.20) 
T po Us 


From the design specifications discussed in the intro- 
Guctory chapter to this thesis, W,ax was selected to be 
70 degrees per second. In order that the output of the sen- 
sor be linear over the range -W..,£W £W),,,,» the value of 
A @ must be constrained to stay within the linear portion 
of the AP versus O curve. Analysis of the cata given Dy 
Sarpkaya has shown that A C) should not exceed roughly +39 
Bor linearity. 6 | 

Inspection of equation (2.20) shows that Ro, Poo? E, 
and Us are still to be determined. R, was chosen on size 
considerations of the sensor and was selected to be three 


inches. The value of Ug will depend on the magnitude of the 


* See Appendix A.l for detailed numerical calculations 
in this section. 


KA 
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desired output signal AP and can be determined from  equa- 
En (2.19). Once Us is determined, « can be calculated from 


the relationship 


y/2 = Lo - 20(2r3) Us (252%) 
where rg > radius of sink 
c - radius of probe 


In order to limit < to a reasqnable value, it is obvious that 
the magnitude of /AP must be kept relatively small. For this 
design a value of AP = 0.1 psi to occur at We Wray was 
selected. 


v 


Having determined y, the xeynolds Number can be deter- 


mined from 
no = 9 no (2.22) 
U R < 
O 
where y = p - kinematic viscosity 


e 


thickness of vortex chamber 


ey 
I 


With the value o? Re determined, the experimental Gata given 

by Sarpkaya on viscous efficiencies can be entered and a value 

of E determined. 
Knowing all the values of the variables in equation 


(2.20) except r,4,, allows one to solve for this variable. The 


po 
position of the probe holes from the center of the sink-tube 
(distance a in Figure 2.3) can then be determined from the 


relationship 








2 2 i 
n a C (anaa 
po EE | 


Merch can be verified by inspection of Figures 2.3 end 2.5. 
Une word of caution must be made concerning a minimum value 
Of Too» Viscosity plays an increasingly important rol» as 


r epproeches zero in a real vortex. If r.o is too large, 


p 
full advantage is not made of the amplifying properties of 
the vortex. On the other hand, if it is made too small, vis- 
cosity rapidly Cecreases the strength of the vortex. This 
phenomena remains to be investigated more fully and is men- 
tioned here only for design consiceraticns. 

The design Cata used in constructing the rate sensor 
Meesumumerized in table 1. The material used in construction 
was plexigless. The porous ring wes mace of felt facec with 
fine wire Screening which gave rigidity to the ring and acted 
as a Slow straightener. “he probe was constructed by cutting 
E 1 3/4 inch piece of number 19 gauge stainless steel hypo- 
Seenic tubing in balf ana inserting & plug. The number 80 
Drill size holes were then drilled 0.095 inches apart. This 
construction allowed the holes to be aligned by rotating one 


side of the probe relative to the other. 


2.5 Signal amplification 


The magnitude of the signal output (pressure and flow) 


by itself is too smail to be used conveniently in the control 
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TABLE I 


SUMMARY OF RATE SENSORI DESIGN DATA 








Specified Data: 


CO max 


Ro 


Â Prax 


AQ pax > 


Arbitrarily Chosen Data: 


Ts 


70° /sec 


o iñ: 


‚0.1 1bf/in* 


30 


0.109682 in. 


0.0212 in. (#419 gauge hypodermic 
tubing OD x ID- 
0.0425 x O02 70 10) 


probe hole diameter = 0.0135 in. (#80 Drill) 


f 
Derived Data: 


Poo = 
a = 
Us = 
Q E 
E = 


0.0512 in. 
0.0475 in. 
171.2 ft/sec 
9.05 ft?/min 


0,50 





system loop and therefore must be amplified. The state of 
the art of pure fluic signal amplification, though not yet 
mature, is advanced to the point where ea combinetion of ex- 
perímental and analytical methods can be used to obtain rea- 
sonable signal emplification with tolerable signal to noise 


ratios. 


The design problems facec at this stege of the develop-. 


“ment of the rate sensor sre: 


1) Given the output pressure-flow reletions of the 
rate sensor probe, whet ere the supply pressure settings 
needed at each emplifier stege to insure linearity of the 
signal? 

2) How many stages of amplification are needed to 


raise the pressure and flow to acceptable levels? 


The lest question could not be answered at this 
stage of the design for the rest of the systems com- 
ponents had not yet been designed. Therefore, it was 
decided to use two stages of amplification and after 
the amplifler stages were attached, to accept the out- 
put pressure-flow relations of the final emplification 


stage as the output cheracteristics of the sensor. 


The first question, although a relatively straight for- 


ward design auestion, is covered only lightly in the literature. 
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The procedure used to obtain the supply pressure settings 
described belcw is by no means exact but it will let the 
Bener Set the pressures close to their optimum settings 
ERCthe first try. 

Pie smpilfvers used to amplify the rate signal output 
were Corning Glass ee Size Proportional Amplifier 
The important dimensions of these devices are: 

power nozzle wicth = W = 0.010 in. 

control rozze wito =z]. 9w 

depch = 2.299 
e receiver distance = 10w 

eure 2.7 is a sketch of the interaction region of a typi- 
Bar proportional amplifier which operates princivly by mo- 
mentun interaction. frollowing the development given by 
Dexter", the deflection of the power Stree to one control 


momentum only is determined as follows: 





V 
tan 6 = E 
Va 
where V. +z, mo Ve. 
Va c Dy M: 
Berk 


and my* mass flow rete of the power jet 
Mo = mess flow rate of the control jet 
t Ve =velocity of the control jet 


V. = Velocity of the power jet 
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Brnical Proportionel Pure lui @ olifier 








SV tanf = -= Me 
P Vp 


but mv = pav? 


4 


| ten = Pac es = | CW Vo" 
2 3 2 
Pa, A Pa Wb Vy a 
where d - depth 
We = control nozzie width 


Wp = power nozzle width 


2 
gc (2.24) 
Wo de | 


N tan O - 


but V^ = (2/9)? 


', ten = We 2/P po = w 


e ? "o Po 
2 
w /P 2. 


(2:595) 
Wo Pp 


Reilly and Moynihant® 


indicate that the output of a 
- proportional device will remain linear provided the deflection 
of the power jet at the receiver ports does not exceed one 


half the nozzle wiéth. Thus: 


tan O = — for linearity 


where L ‘= receiver distance 


Combining this with equation (2.25) gives: 


w > He Pe for linearity 
, 2L © Wp Pp 
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$ 
4 , 


INS result iS readily extended to two control ports with the 


result: 
_"p_ > Wo lrorPer) for (2.26) 
eL Wp Py Linearity 


Having devéloped the basic equations of the amplifier, 
the output characteristics of the rete sensor probe must now 
be investigated. According to potential flow theory, if the 
probe is set so that the probe holes are at 45 degree angles to 
the oncoming flow, the pressure in the probe chambers is. below 
ambient. Since the probe chambers are connectec to the control 
ports of the first stage amplifier, this means that only the 
pressure differential existing in the two chambers anà com- 
municated to the first stuge control ports is available for 
deflecting the jet. This is undesirable as it results in 
less gain when compared to operation with positive flow through 
both control ports. nowever, due to viscosity the ideal flow 
pattern is radically altered so that the pressure existing Jon 
the chambers when the flow angle is set to 45 degrees is posi- 
tive. 

The various flow rates to be expected out of either 
end of the probe, as the main sink flow is deflectec, can be 
approximated by setting the sink flow at its operating value 
anc turning the probe with respect to the flow. The results 
of such an experiment are shown in figure 2.8. 

The output pressure-flow characteristics for a given 
Py are shown in Figure 2.9. An infinite number of such curves 


exist since P, varies with A Q- 
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Probe 3oundary Pressure versus io Load Flow 
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figure 2.9 


Probe Onutnut Prosare- “low Characteristics 
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The procedure initially used to set the first stage 


Supply pressure (P L was to Getermnine the output pressure- 


Y 


+ 


Slow relationship for maximum AP and superimpose the input 
characteristics of the first stage amplifier to determine 
= Az, 


supply pressure was determined. However, when this procedure 


Then, from equation (2.25), the 


>> >) : 
Mer y Pon) max nax* 


was followed, it wes found that, because the contrcl port area 
Of the first stage anplifier was so lorge compared to the probe 
hole area (nearly three na large), the input impedance 
of the emplifier wes nearly zero. 

Circumventing this problem, equation (2.24) was used 
in Memeor (2.25). With. the constraint placed on tan @ and 


for two control port operation, equation (2.24) becomes: 


Wo > Wo (vane = ee (2.27) 
SL Wp Vo? 


From measurements of the flow rate occurring at the 


expected maxirum pressure just outside the probe hole, 


2 5 i i e x e = 
E" - Vor ee was obtained fron VeL Van = 


C 
to be obtained fron 


(V 


(QoL mex T LeR min) Which then allowed V 


equation (2.27) 


P 


To obtain Pp the following relationship was used: 


u 


è + Cap Ap VE E, (2.28) 


p dp 
.where Lo - Ap Vo 
C discharge coeficient of the power 


dp nozzle which was determinec from the 
design data sheet published by Corning 
glass works. 
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was determined to 


4 


From the calculations in Appendix ÀA.2, Fo 
me 7.7 inghes of water. 

The second stage pressure was set using informetion 
given in the design data sheet publisheà by Corning Glass 
works for these amplifiers. This sheet states that for max- 
inum gain the control bias level pressure shoulc be one tenth 
the power jet pressure*. 

The bias level of the second stege controls wes deter- 
mined by setting the first stage supply pressure at 7.7 
inches of water anc measuring the pressure in the line connect- 
ing the output of the first stage anc the control port of the 


first stage. This procedure resulted in an initial setting of 


Poe equal to 5 psi. 


2.4 Rate Sensor Dynamics 


This rate sensor is cheracterizeä Are 1) 
a pure time delay resulting fron Ehe time it takes for the flu- 
id particles containing tangential momentun to reach the probe, 
2) a first order lag due to the resistances of the probe holes 
and control nozzles of the first stage amplifier and the vol- 
ume connecting them, and 3) the dynamics of the staging am- - 
plifiers and the lines used in interconnection. 

The transport delay is the most significant. It may 

$ 


* This procedure used in setting the first stage 
pressure resulted in a Pol = 1.95 inches of Water and hence 


was not used. 
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p" 


De approximsted by the time it takes for the sensor to con- 


pletely renew the supply of fluid in the vortex chamber 


ane sink-tubes. Thus: 


TEs 


* 
= 


where V 


Q 


The totel volume (V) 


y z 


2 


mL seconds (2.29) 


Q 


total volume 


flow rate given by eguation (2.23) 


is given by: 


avRofh + 24Yrglyg 


where Ro, h and rs are as before 


qo see 


distance to probe measured along the 
axis of the sink-tude 


For the geometry and flow rates of this device, the value of 


T is 0.0556 seconds. 


The time constent of the first order lag referred to 


in item 2) was found to be extremely small anc was neglected 


in arriving at the transfer function of the rate sensor. 


the dynamics of the staging amplifiers are discussed by 


Belsterling and Tsui® and Boothe’. 


These works show the ce- 


Vices to respond very much faster than e time delay of this 


rate sensor. In addition, the design data sheet for these 


amplifiers indicate thet the frequency gain is over 1000 cps, 


well out o? the range'of the rate sensor. 


Thus the transfer function relating the output (AP) 


to the input anguler rete (W) is given by: 
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pu Ae 
E l 
Where - Ay = steady state gain 
T4 - pure time delay given by 
t equation (29) 
3 = Leplace operator d/ét 


EE Testing Frocedure onc results 


the rate sensor was made into a self contained pack- 
ege by atteching the staging amplifiers directly to one side 
of the vortex chember. The sensor was equipped with a man- 
ifold allowing air to be supplied with only one line. This 
E aroa pressure was then reduced to the required pressures 
needed at the vortex chamber supply plenum and power jets 
of the staging amplifiers. Variable resistors were used to 
set the pressures and were constructed from tapered cowel 
pins through which holes were drilled.  PLhotogrsphs of the 
assembled rate sensor are shown in figures 7.4anä D.l . A 
differential pressure to electric transducer was attached to 
the outputs of the second stege amplifier anā the assembled 
đevice wgs tested on a rete table loi TUE 
Laboratory of the Massachusetts Institute of Technology. 

The testing procedure consisted of the following: 

Y) Steady state test. over the rangeul = nio 
racgians per second with the output loaded by orifices approx- 
imating the size of the control parts of an amplifier similiar 


to those used on the rate senscr for signal amplification. 





| í 
2) Steady state-test over the seme range with the out- 


put flow equal to zero; i.e., infinite load. 
5) Frequency response test at discrete frequencies 
up to 4.67 cps. 
4) oteady state test while the device was being simu- 
taneously rotated about its sensing axis and an arbitrary axis. 
Ihe test results for the first three tests are shown in 


Figures 2.10, 2.11, 2.12, and 2.13. 


2.6 Discussion of results 

The results of the steady stete rate tests show the 
following: | 

l. The very small signals available from the probe can 
em@ecd be amplified into reasonaodle signals. 

2. ihe output of the sensor is not quite linear over > 


Sted. 


xD 


the range t 
3. The gain of the sensor with ciockwise rotation is 
Significantly larger than the gain with counter clockwise ro-: 
tation when the bias is subtracted out of the results. 
4. The gain incre&éses continually as the input rate 
is increased. : 
5. The dog leg effect produces a highly undesirable 
rapid change of gain in the counter clockwise direction. 
In an attempt to'explain the last four items, accurate 
flow and pressure measurements were taken on each side of 
the probe as the angle of incidence of the probe holes with 


the oncoming flow was varied. Uriginally, in aligning the 
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Figure 2.11 


Rate Sensor Steady State Output 
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Typical Vortex Rate Sensor Frequency Response Data 
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t, 
probe holes (the turning of one probe hole with the other helå 
fixed), only the dead-ended pressures were used, taking the 
probe holes as being aligned when the pressures were equal. 

Às was pointed out in article 2.3, signal amplifica- 
tion is essentially dependent on the flows coming from the 
probes. ‘shen these flows were measured, it was found that 
equal dead-ended probe pressures resulted in radically aif- 
ferent flows indicating that, although the holes were sup- 
posedly of the same area, they were behaving as if they were 
not. 

The probe was removed from the sensor and observed under 
a 50 power microscope. This showed that during the machining 
process, the drill had been allowed to cut a trough into the 
entrance to one of the holes. The effect of this trough was 
apparently to catch more flow from the main stream resulting 
in different orifice characteristics. 

The effect of having the two probe holes producing 
different flows for equal boundary pressures was as follows: 
It has previously been shown that, as the sink-tube axial 
velocity is deflected, the pressure outside the probe holes 
increases or decreases depending on the sense of angular ro- 
tation ana the probe hole in question. As the pressure out- 
Side the probe hole which gives less flow for a given pressure 
is reduced, a point is reached where the flow approaches zero. 
By referring to Figure 2.8, it can be seen that, as the flow 
approaches zero, the relation between the pressure tes 


the hole and the resulting flow becomes highly non-linear. 








This highly non-linear low region prouces the dogleg ef- 
fect noted above. This effect is not present when the direc- 
tion of rotation is reversed for, as the boundary pressure 
outside the hole which gives the greater flow for a given 

Fi Sure is reduced, the non-linear portion of the curve is. 
avoided*. At the same time, the other probe hole íis experi- 
encing inoreasing pressure and hence this combination prođuces 
no dogleg effect. Thus, to remove the efrect, the probe hole 
without the trough was altered to give approximetely the same 
pressure-flow characteristics as the other probe hole. 

To explain the difference in the maximum signal output 
between clockwise and counter clockwise rotation, one has only 
to realize that the positioning of the probe holes with respect 
to the center line of the sink-tube is critical. If the holes 
are off center the slightest amount, the output will be de- 
creased on one side and increased on the other resulting in 
different peak outputs. 

The effect of the gain continually increasing with in- 
creasing Wcan be partly explained by realizing that deflec- 
tion of the first stage power jet is roughly proportional to 
(A). “Generally the output pressure of a proportional am- 


plifier is given by 
AND. 


KAP, 


Dae ren K, (A <) 





* At least over the range of rotations tested. 


lu 





hence AP, = Ks VA 


Now as shown in Figure 2.8, Aw is approximately proportional 
to W over a substantial renge of operation, making the output 
pressure generally Du 5 oe to the square of @. This 

Of course is only approximate but does roughly explain the 
increasing gain of. the device. In fact, if the output pres- 
sure is converted to flow using the output characteristics 

of the Canna se 3 one finds the output flow quite 
linear with respect to QJ. 

Although accurate verifications of the explanations | 
given for items 2, 3, 4, and Sh bore were not accomplished 
due to inaccessibility of the rate table, the explanations 
given were partially verified and showed that the dogleg had 
been removed but that Se improvement was made in 
balancing the outputs by centering the "o 

The zero input bias signal was found to be oniy partly 

«Gue to the unequal flows from the probe. The amplifiers 
themselves are biased when AS control flows ere zero. This 


can be reduced somewhat by proper arrangement of the cascade 


so that the output bias of one tends to cancel the output 





bias of the next stage. However, adjusting orifices general- 
ly must be placed somewhere in the cascade to remove the bias. 
The results of the frequency response tests verified the 
| pure time delay nature of the dynamics of this sensor, The 
experimental phase data is approximate at the low frequencies 


because the method used to obtain the data (the input and 
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output were recorded on permanent recorcing paper and were 
compared in order to measure phase) coulé not discern small 
phase shifts. The amplitude data, however, is quite accurete. 
Lissajous figures on an osciloscope could not be used 
to obtain frequency response’ data because the output noise 
moved the figures sround the scope precluding any measurements*. 
The best data was obtained by taking pictures of the 
oscilloscope traces and measuring phase ana attentuation data 
from these. The rapid increase in phase verifies the existence 
of a pure time delay. The deviation from theory in phase 
shift and attenuation is probably due to volumes in the lines 
connecting the first and second stage amplifiers and the vol- 
ume and resistance of the orifices between the output of the 
second stage amplifier and the pressure transducer. 
The results o? the last test showed the et signal 
did not change while the sensor wss being simultaneously 


ERS 
rotated about its sensing axis and an arbitrary axis. 


2.7 Conclusions and recommendations 

Throughout the cevelopment of this rate sensor little 
attempt has been made at optimization. The objective was to 
develop a rate sensor capable of being used in a control systen. 


To this end the following conclusions are drawn: 


* The not'se of fluidic proportional amplifiers is of 
a peculiar nature. Although high frequency noise iS prevalent 
as shown in Figure 2.13, a more cisturbing very low frequency 
noise is present which causes the signal to wancer and make 
unexpected excursions which appear almost as steady state bias 
changes. 
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I. USG Or tne design procedure suggested in article 252 
will allow the designer to theoretically Gesign a rate sensor 
peculiar to his applications. 

2. The steady state gain of the unamplified signal 
may be predicted analytically but due to the low input in- 
pedance of the type of amplifier used, the amplified gain will 
have to be determined experimentally. 

3. The dynamics of the sensor may be determined to a 
reasonable degree or accuracy by theory. 

4. The output of the sensor is very dependent on the 
accuracy of the Arne of the probe holes and alignment in 
the sink-tube. 

5. Noise of the very low frequency type is a serious 
problem causing unexpected wandering of the DC level of the 
signals. 

6. Contrary to the results of sServkaya, it was found 
that closing one of the sink-tubdes and keeping the flow con- 
stant, resulted in an increased output signal for a given QJ. 
Since this eouivalently means nearly the same signal with less 
power consumption, the author recommends thet sensors of this 
type be built with oniy one sink-tube and, as recommended by 
Sarpkeya, that the end opposite the oed ae be designed to 


streamline the flow into the sink-tube. 





CHAPTER III 


ANGULAR POSITION SENSOR ANALYSIS AND DESIGN 


f 


l Criteria Used in Selecting Operating Principle 


Angular position of the system can be obtained by 
ane log® or &igital?»? integration of the angular rate signal 
available from the vortex rete sensor. Digital integration 
woulda possibly involve developing a pressure sensitive os- 
cillator and some form of digital-to-analog converter, both 
of which sre beyond the time limitations of this thesis. 
Since a linear position signal over a large range of angular 
positions for demonstration purposes was cesirec, analog 
integration was ruled out since very lerge ‘volumes would have 
been required to"store"the position sivhal. Furthermore, the 


"7 


State of the art of pure fluid analog computation’ is such 


that up-down integration is inaccurate. : 


4 
y 


Eoo Position sensor Description 


Because of the problems mentioned above, it was decided 
to use a moving part position sensor oonsisting of a penáulum 
set between a pair of flapper nozzle velves to obtain a dif- 
ferential pressure signal that is a function of angular posi- 
tion. A schematic of this sensor is shown in Figure 5.1 and 


a photograph of the final design is shown in Figure D.d. The 
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. igure 3.1 


owetcl cs the inculber +ositiox Sensor 











device is designed so that as the sensor rotates, the pendulum 
mrys to Return to a vertical position. In so doing, it widens 
the gap on one side anc closes the gap on the other until the 
moments caused by the flow forces from the nozzles, balance 
the mcment due to the weight of the pendulum. «Jhen this 
occurs a differential pressure exists in the chambers and 
constitutes the output signal of tne sensor. 

It is readily seen that the output differential pres- 
sure signal varies as the sine of the angular position and 


that the sensor also senses direction of angular position. 


Meo Discussion of the Dynamics of a Pendulum Angular 


Position Sensor 


If the bearing which supports the pendulum is relative- 
ly free of damping, the angular position of the pendulum is 
unstable and the slightest perturbation will start the device 
oscillating. in the bang-oany control envisioned for the 
final system, this oscillation is undesirable. To understand 
the reasons for this instability, a lineer model of the system 
can be constructed to describe the behavior of m system 
when the angular position of the pendulum is smell. Figure 3.2 
shows an anklytical model of the sensor. Summing moments 


about the pivot point gives the following differential equa- ' 


wion: 


2 e ° 
E 1,02, Ot (01, + WA = Tato) 





Figure 3,2 
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e 
Ps 


o 


where 0 ag/at 
Dl» equivalent viscous camping 


K e;uivalent spring constant of the 
flapper nozzle valves 


Ta a disturbance torque 


Neglecting the term WO, since it is much smaller than 


xlna O» anê introducing the symbols J= > Ne. B= bl,, 
K= klp, and s=-d/dt gives 


E + ES + 3 O = jt (5.1) 
The spring force K Mis given by 
K O (P,à - Pgàllh —* Aral, 
where P, = pressure in one chamber 
Po pressure in the opposite chamber 


rearranging gives 





e. (3.2) 
The linearized transfer function between AP and @ 


is, in general, a first order lag of the form?®, 


Al 


eS 
" 


n 





2n (3.3) 
6 Ts+l l 
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Substituting equation (3.3) into (3.2) gives 


| . 
x = Ala = C2 (324) 2 
T 391 Ts+l 


Ee this equatien is substituted into equstion (3.1) the 


result is 





ES + Bs + 92 | 9 = Ti 


Tsrl 


From this the transfer function between 9 and Ta is 


E z Ts +1 
! v. 
Ta IT s°+(IJ+BT)S2+ BS&Co 


The characteristic equation of this transfer function is 
Lu + (X4 BT] s^ + ES 4 3 


The Routh stability criterion for linear systems gives the 


condition for stability as 


B(J4BT) > JTC 


Or 
Pa se l 
B> N Ira Gp a (3.5) 
27 


where the positive sign for the radical has been chosen since 


the damping can not be negative. If equation (3.5) is re- 


B> 1/2 J|“ £ 4JC2 = g 
m T 


E 


written as 


e 


ts 


^. 
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it can be seen that the value of C> ís the important para- 
meter determining the stability of the system. As Co in- 
creases, B must be made corresponcingly larger. From 
equation (3.4), Co is the value of the steady state spring 
constant of the flapper-nozzle valves. Thus, as this spring 
constant increases, the magnitude of the damping must be 
meer eased. 

Initially, to introduce damping into the system, the 
ball bearing used to support the pendulum was packed with 
heavy gresse. This did not introduce a sufficient amount 
of damping into the system. The problem was finally sa 
by placing two small squeeze film dampers on either side of 
the pendulum end adding a minute amount of grease to aid in 


A 
the damping process . 


9.4 The Influence of Acceleration on the Output Signal 


The ‘output signal of this sensor will be influenced 
by the angular acceleration o? the device to which it is 
Ë achod. Figure 3.5 is a sketch showing the position sen- 
sor mounted on the system and offset from the axis of rotation 
by a distance r. 


if the sensor is given an angular acceleration q , 


* Over most pressure settings the air in the squeeze 
film dampers provided sufficient damping. However, if the 
chamber pressures were made too high then the size dampers 
used were inadequate. 
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then the end of the pendulum receives a linear acceleration 
Po. The torque on the pendulum is given by 


* 
Ta fe e - ra (3.6) 
& 


This torque must be balanced by the pressure forces on the 


enc of the pendulum giving 
DES 4 ra 
Pr 


ND = wr (Gal) 
LA CL 
gal, 





This AP corresponds to an error signal. It may be reduced 
by making r small assuming the other parameters are fixed. 


“The moment of inertia (1) of the pendulum is given by 


Thus the resulting torque on the pendulum is given by 


T = —- iC (3.8) 


This torque too must be balanced by the pressure forces on the 


end of the pendulum giving 





r - Vi 
A Pa Aly = — Y, G 


or 


2 
Pye = A a (3.9): 


* ‘Subscripts LA 6nd ña refer to linear acceleration 
and angular acceleration respectively. 
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This Â P also represents an error in measuring angular position. 
The total error signal is the sum of the two differential pres- 


sures given by equations (3.7) and 3.9) 


, e 2 
AP AE. + Tala (3.10) 
error Y gAl, | Axe ; 


This pressure error may be converted into an equivalent angu- 


lar position error by the equation 


0 = sin? E Al, (G 13) 
Wy 
which is obtained by a steady E torque balance on the 
pendulum when at en angle@ from the vertical. Substitution 
of the value of AP from equation (3.10) into equation (3.11) 
Will give the error in terms of anguler position. Primarily 
because the accelerations to be ceadt with in this system are 
quite small, substitution of the actual values into equations 
(3.10) and (3.11) showed that the maximum probable error would 
be less then one degree. if this system were to be made into 
a very fast high torque system, this type of position sensor 
woulé probably be undesirable unless the decer vishbd to 


have acceleration feed back in the system. 


3.5 Transfer Function 

Because the amount of damping added to the system to 
make it stable is difficult.to determine, it was decided to 
connect the output signal into a pressure to electric trans- 
Gucer anc experimentally determine the transfer function. By 


e 


Observing step responses on an oscilloscope, the system dynamics 








can readily be determined. 

The procedure was as follows: 

1) The eu me signal was connected to the transducer 
and loeded into orifices which approximated the size to be 
expected at the control ports of the summing device. 

2) The output signal was then displayed on an oscil- 
Boscone. 

3) The pendulum was given an initial deflection and 
then released. 

The response of the pendulum was observed to be es- 
Sentially that of a first order leg. By measuring the time it 
book the signal to decay to 36.8 percent of its initial value, 
gave a close approximation to the time constant of the sensor. 
By averaging & large number of such measurements the time con- 
Stant of the system was found to be 0.04 seconds. This time 
constant was found to be much larger than the time constant 
associated with the flapper-nozzle chambers". Apparently, 
the measured time constent is essentially that of the mass- 
damping time constent of the systen. | 

A typical steady state output curve is shown in 
Figure 3.4 from which the steady state gain may be measured, 

The: trensfer function, relating the output pressure 
Signal to the input angular position signal for small valves. 


of 6 , is given by 


* The time constant cue to the volume of the flapper 
nozzle chambers and connecting lines was found to be 0.008 
Seconds. The calculations are not shown here as they can be 
found- in any standard reference such as reference 10. 
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where Ko = gain to be determined from a plot 


similar to Figure 3.4 once the final 
system pressures have been set 


T2 0.04 seconds 


Before concluding this chapter it is to be noticed 
that little information has been given concerning design 
procedures. This was purposely done as such procedures are 
more acequately covered in the litereturel9, However, the 
designer must keep in mind the over all system ano the re- 
lationship of each component to the others. For example, 
it is obvious that the output of the positicn signal is to 
eventually be compared with the rate signal. It is therefore 
reasonable to assume that these signals should be of the 


Seme order of magnitude. 


bro switching Curve Computer 


As shown in Appendix C.1, the ideal switching curve 
for this system is a parabola passing through the origin. 
Bue slope of this curve is infinite at the origin. Since it 
was planned to generate the switching curve by operating on 
the output signal of the position sensor, this meant thet as 
ĝ varied slightly from its zero position an infinite gain 
mOuld be required. Because of this consideration and the 
fact that the output signal of the rate sensor is not depend- 


able at very low angular rates, the output of the position 


hy 


D 





J 


or, without modification, was used as the switching curve. 
Such a sine curve output can be adjusted to approximate 

| parabola over a large range of position signals as shown in 

Figure 3.4. The approximation is poor when 6 is small, how- 


ever, and thus, it is to be expected that "linear switching" 


will take place when the errors introduced into the system 











CHAPTER IV 


SIGNAL COMPARATOR 


4.1 General 

The method! selected to algebraically sum the pressure 
signals from the rate sensor and position sensor consists of 
adding additional sets of controls to & proportional amplifier 
end having the power jet deflect in accordance with the vector 
. sum 0f their momentums*. A stream interaction proportional 
emplifier, Operating on this principle, is available from 
Corning Glass Works** anà has two sets of control ports. The 
design data sheet accompanying this amplifier does not in- 
dicate that the second set of controls can be used for sum- 
ming and, therefore, gives no data on the gain of these ports. 
In fact, the ports sre not equipped with fittings and, there- 
fore, had to be adapted for use in this control system. A 
photograph of this "summer-amplifier" is shown in Figure D.Z2. 
The vertical controls are the seconcary controls and were’ 
adapted for use by inserting hypodermic tubing into the holes 


and then cementing the tubing into place. 


, 


* Other methods of summing signels, using no moving 
parts, are available and could have been used in this system. 


** Catalog Number FL2511-2-1211 September, 1965 
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ENS Exverimentel Analysis 


À silhouette of the emplifier described above is given 
mo the cesign data shect. From tnis Silhouette, it was found 
meat the width of the seconcary control ports is less than 
the wicth of the primary ports. In eddition, the secondery 
ports are not perpendicular to the power jet. Bi oo these 
geometrical difrerences, it is obvious that if the same pres- 
Bare differential is applied to both sets of controls ports, 
the power jet will be deflectea a greater amount by the pri- 
mary control ports. As a consequence of this greater deflec- 
tion, the output pressure will be larger, Therefore, the 
primary m ports will exhibit a greater gain. 

To determine, the difference in the gain of these two 
Sets of control nozzles, a number of tests wore performed on 
the amplifier. The tests consisted of l) systematic varia- 
tion of the control port bias levels about the levels ex- 

pected from the outputs of the rete and position sensors, 

2) systematic variation of the supply pressure, and 3) at any 
particular bias level and supply pressure setting, recording 
the gains of the two sets of control ports. 

Typical data from these tests for one particular com- 
bination of supply pressure and bias level settings is shown 
in Figure 4.1 During all these tests, the output of the 


summer-amplifier was losded into a manometer. 


4.3 Discussion of Results and Conclusions 
The results of these tests showed that the gain ratio 


of the two pairs of control ports did not remain constant as 
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the control bias level or as the Supply pressure was changed, 
However, for any given set Of control bias pressures and a 
given supply pressure, the ratio of the gains of the two con- 
trol ports could repeatedly be obtained after signals were 
introduced and then removed. 


The results further showed that, once the ratio of gains 


had been determined at a part iculan- setting of bias levels and 


supply pressure that the algeoraic sum of two control signals 
could accurately be predicted. | 

From the results of these tests, it was concluded that 
this device could be used as a Summing device with accurate 
and predictable results if the bias levels df the input sig- 


nals and the power ‘jet pressure were known. . 











CHAPTER V 


4 
DRIVING BISTABLZ AMPLiIPlixs w1lTh PROPORTIONAL AMPLIPLERS 
f 


5.1 Introduction to, Input Control Impedances of Bistable 


Amplifiers?® 
in Figure 5.1 a typical bistable fluid jet amplifier is 


shown with the power jet attached to the left boundary wall. 
Mae this condition, the entrainment flow through the left con- 
trol port is greater than that through the right when both 
control ports are open to the atmosphere. if the left con- 
trol port is sealed from the SUE while the right port 
1S left open and a plot OT the lert control pressure and 
flow is made, a plot similer to Figure 5.2 will be obtained. 
This graph shows that as the left control pressure is in- 
creased above atmospheric pressure, a point is reached where 
the jet is supplied all the entrainment flow it can take. It 
then becomes Separated from the left wall ana switches to the 
right. hen this occurs, the flow sudcenly decreases Que to 
a different pressure in the interaction chamber and now pre- 
sented to the left control port. vil the lett control por 
pressure is then decressed, a point is reached where the jet 
will switch back to the left side. | 

If the amplifier has perfectly symetrical control ports, 
the "stream-on-right" curve will represent the pressure-f low 


characteristics at the right control port when the stream is 
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attached to the left wall. From the graph it is seen that the 
flow differential required for switching from left to right is 
Qa7Qq and the pressure differential required is Inn 
EN larly the flow differential to switch from right to left 
is Qp-Q,, and ene qus pressure differential is Patm 3 D. 
The two sets of flow and pressure differentisls are not nec- 
cessarily equal. 

The effect of loading on the output legs of these de- 


vices is generally to decrease the flow and pressures required 


for switching. 


5.2 Problem Statement end Objective of Experiments 

Because of the nature of operation of proportional 
amplifiers, the useful signals from these devices are elevated 
about e bias pressure level greater than atmospheric. © It wes 
mot clear what the effect of this bias level would be on the 
switching characteristics of a bistable amplirfier being driven 
by a proportional amplifier. This information was necessary 
before undertaking the design of tho bistable reaction Jets 
in order that the pressure differential required at the con- 
trols of the proportiòonal amplifier, to effect switching of 
the reaction jets, might be minimized. 

The test procedure vas to bias the control ports of a 
memory device” at various pressure levels above atmospheric 


anà then measure the pressures anà flows recuired to effect 


switching at these elevated levels. 





^ A memory device was selected for these tests as this 
is the type of device used for the first stege of the reaction 
jets discussed in the next chepter. 
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During the tests the output legg- were loaded into an 
rifice whose area was approximately 0.75 times the receiver 


areas of the device. 


Meo Test Results and Conclusions 

The memory device selected for these tests was a 
Corning Standard Size Load insensitive EL 

A graph of representative results is shown in Figure 5.93. 
It can be seen that the effect of the bias level is to radi- 
cally increase the flow and pressure differential required for. 
switching. f 

During these tests it was found that the input charac- 
teristics of Bach control port were quite different.  There- 
fore, to more accurately show the flow and pressure Giffer- 
entials required for switching, the pressure-flow relations 
for both legs M plotted. 

when a bias level is present at the control ports of 
a bistable amplifier, the power jet is being supplied more 
flow than it can entrain. In this condition the bias flow 
is attempting to operate the device in a proportional manner. 
However, due to the close proximity of the bounding walls, 
the jet will still reduce the pressure slightly on cone side 
anc attach itself to that wall. This increased pressure ° 


differential required for switching can not be tolerated in 


this system as it represents a large hysteresis in the 





* Catalog No. Fb 2212-2-1211 
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switching of the bistable reaction jets. This problem, how- 


eyer, can be eliminated by making the bistable amplifier 
large enough so thet the entrainment flow of the power jet 


is just equal to the bias level flow. 








CHAPTER VI 


Semon AMPLIFIER aND BISTA9LE REACTION JETS ANALYSIS AND DESIGN 


65.1 Determination of the Required Output Thrust 


So far little mention hes been made of the dynamics of 


the system to be controlled. As mentioned in the introductory 





chapter, the system was to approximate an ideal inertia load. 


The dynamics of such a system are given by 


È 


10 (6.1) 


where T “= the applied torque to the system 


ey 
I 


polar moment of inertia 


D 
I 


angulsr displacement 


differentiation with respect to time 


If this equation is integrated once with respect to 
time and the initial angular rate is assumed to be zero, one 


obtains 
03: was Lx (6.2) 


Integrating equation (6.2) with respect to tíme and setting 


the initial angulsr displacement to zero gives 


2 
Ó - = D | (6.3) 


L > 





rom the design considerations discussed in the in- 
0 
nr QC tO wW Chapter, it was cesirecd that the systen respond 


a 


De able GO receive large ansuler aissiuce= 


rt 


Ey siowly and 
ments. On the basis of these consicerations, it was cecided 
that the parameter T/J shoulc be such that if the system 
Starts from rest under full acceleration, Ít must be travel- 
ing at an angular velocity equal to 1.1 radians per second 
et the ena of 45 degrees of angular displacement. 

If equation (6.3) is solved for t ané substituted 


into equation (6.2) the result upon rearranging is 


HN TIE SCE (6.4) 
g^ 20 


Substituting the values of (O aná C) into thís equation gives 
ET - 0.772 sec”. From equation (6.3) the time to travel 
5°, with this T/J ratio, is 1.424 seconds. 

Having determined the T/J ratio, an estimate of the 
systems finsi polar moment of iE vias made oy accurately 
determining the polar moment of inertia of the rate sensor and 
Eultiplying the result by a factor of three to account for the 
remaining equipment.to be added to the system. Tne J so cal- 
CGulatec was foung to be 7.5 in-ozf-sec^ "wniechomesults uua 
required torque of 5.79 in-ozf. lf the lever arm is assumed 


to be 5.79 inches, then the thrust requirec is one ozf. 


) 


6.2 Introductory, Theorv of staging Unvented Bistable 


Amplifiers? 13, 14 


This system was to Gevelo» the required thrust by 


using the change of momentum in a renction jet. since the 





i 





# 


power available from the output legs of the summer-amplifier 
Was too Smell to develcp the required thrust, this signal 
mec. tO be amplified. | 

Compagnuolo*” has shown thet this power amplification 
can be efficiently accomplished by staging bistable emplifiers 
in ascending order of size with each stage operating at the 
same Supply pressure. He has pointed out that in orcer to 
meke such a cascade, two design principles must be adhered to: 

l) For maximum pressure recovery, the splitter ĉis- 
tances must be kept small, preferably in the core region of 
the power jets which in general is approximately six nozzle 
widths down stream. | 

| 2) E order oe the units be capable of being loaded, 
each stage must have a high degree of stability which implies 
moving the splitter Gown stream to over approximately twelve 
nozzle widths. 

At the present state of the mU a bistable unit can 
not possess both the cherscteristics of high pressure re- 
covery (high efficiency) and high stability. As the splitter 
is moved down stream, the pressure which can be recovered, 
falls Nes if the splitter is moved up stream, then the 


ability to operete under load, (stability) falls off rapidly. 


" Tne core region of the power jet is that region in 
which the velocity of the jet has not been reduced from its 
value at the power nozzle due to viscous friction with the 


Surrounding fluid. 
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Compagnuolo suggests that since the last stage is the 
most important for maximum efficiency that it be mace with the 
splitter located in the core region of the power-jet. This 
stage can then be held into the load by making the preceeding 
Stages high stability devices. The result is that the overall 
cascade exhibits high efficiency and stability under load. 

In general the instantaneous flow gain of high stability 
devices is approximately ten. This means that it is possible. 
meo increase the size of the power nozzle of each stage by 
roughly a factor of ten. lit should then be possible to have 
flow gains in a three stage device of at least 1000 and, if 
the efficiency of the cascade is good, power gains of greater 


than 5000. , 


6.3 Initial Experimental Design 


Figure 6.1 shows an output leg of a bistable amplifier 
turned through a 90 degree angle. If a control volume ís 
drawn as shown and the momentum equation for steecy state 
Operation is written for the x-direction, the resulting 


thrust cue to the exiting jet is found to be: 
e re 
T = PAeVe (S79) 


ASSuming uniform velocity over the exit area, equation (6.5) 


can be rewritten as 


T =P dee | (6.6) 
Ag 


from this last equation, knowing the thrust desired and 


qe 





control volume. 





Pigure 6.1 


Control Volume for Determining the Output Thrust 











assuning a reasonable value for the exit area, one can Geter- 
Mine the @low required to develop the desired amcunt of tnrust. 
If one further knows the differential flow required 
for switching the first stage, then the flow gain needed can 
be readily determined. 
it was shown in the preceecing chepter that it is ce- 
Sirable to have the entrainment flow of the first stage 
approximetely equal to the output bias flow of the driving 
Proportional amplifier. ‘This entrainment flow depends on 


12 It 


the geometry of the amplifier and the supply pressure. 
can be seen thet if sufficient experimental Cats were avsil- 
able that it wouio be possible to 

1) determine the supply pressure required for a chosen 
first stage amplifier design from the entrainment flow versus 
supply aan: characteristics and 

2) with the supply pressure chosen aS ebove, deter- 
mine the differential flow required for switching. 
It would then be a simple matter to determine the flow gain 
requirec and the number of amplification steges needed. 

Lacking any such quantitative Geta, a preliminary 
two stage Seren amplifier wes constructed in order to 
obtain the following information: (listed also are the rea- 
sons ‘this information was desired) 

1) Entre inment flows of a first stage design at vari- 
ous supply pressures to determine the supply pressure setting 
to match the output flows of the summer-ampiifier. 


2) Stability under losd of the unit to check the set- 
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DP. 
5 


ting of the splitter distances, wall set back distences, etc. 
3) Output pressure-flow relationships i determine 
the optimum loading for meximum torque output. 
4) The number of stages required to develop the 
desired torque. 
5) The thrust output for various supply pressures end 
loeds. 
The important dimensions of the amplifiers used 
š 


in this initial design ere as follows:* 


é 


- First Stage 


Fower nozzle Width = w = 0.020 in. 
Control noz2i€ wrath s w 
wall set ba ck = SW 


Splitter distance = 14w 


1 


Aspect ratio = h/w = 1.5w 


Receiver width = Jw 
second Stage 


Power nozzle width = w = 0.0625 in. 


Control nozzle wicth w 


Wall setoack = l.4W 


* Many of the dimensions given here, auch as wall set- 


back and receiver width, are based on carefully considered 
data given in references 12, 12, and l4. Time ana space lim- 
itetions do not allow a full discussion of the effects of 
varying these parameters. 


25 








Splitter distance 3 6w 
aspect ratio , = 1. 999% 


Keceiver width = J ¿OW 


i 


The devices were made out of piexiglass using a pan- 
tograph engraving m-dhine. 

The results of tests performed on this initial design 
showed thet at a suppiy pressure of two psi the entrainment 
flow in the attached side was approximately 0.04 SCFM. This 
flow was approximately equal to the no losá output flow of 
the summer-amplifier when operating at the supply pressure 
expected for final design setting. 


This device was quite steble provided the output load 


y 
4- è 


was relatively small. Typical output curves for the device 
are shown in figures 6.2 end 6.3. in rigure 6.2 the device 
is not being criven by an external signal anā therefore can- 
not be heavily londe. The overload area represents condi- 
tions under which the second stege power jet commences to oscil- 
late. In Figure 6.3. the device is being Griven by another 
bistable amplifier. It is apparent that the device can now 
be more heavily losded before the overload region is reached. 

In order to determine the value of the exit nozzle 
area for maximum thrust, a number of nozzle load curves were 
drawn on a plot similar to Figure 6.5. It was found that 
maximum thrust output of one nozzle corresponded closely to 
the point of maximum power transfer between the amplifier 


cascade and the nozzles. since loading decreases the flow 
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out of one leg, part of the flow is diverted into the other 
eg where@it develops thrust in the Opposite direction. There- 
more, considering the total thrust ceveloped frcm the device, 
it was found thet the point of operation for maximum total 
thrust ley slightly up the output curve toward the region of 
higher flows. | 

The results also showed that in order for this two 
stage device to develop an ounce of thrust that the supply 
Bressure would have to be over ten psi. Since this was ín- 
@emoavible with the entrainment flow desired, the rinal 


Gesign would have to have at least three stages. 


6.4 Final Three Stage Amplifier Tesien 


In order to use the «cata gainec from the two stage 
amplifier to full edvantage, it was decided to use these 
Moret two stages in Se final design and alter them only 
slightly in order to acccmmocete the thira stage. This 
further alloved the power nozzie area of the third stage to 
be calculated from the data obtained on the two stage device, 

The procecure used to determine the dimensions of the 
third stage are as follows: assuming incompressible flow”, 
the flow from the thirG stage power nozzle is given by 


Qs z Ca A 33 Cona Ps (6.7) 


p 


* This assumption is reasonably valida as the final 
device was to Operate at two psi. 
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where Q;3 = flow from third stage power nozzle 
è 


A 43 area of third stage power nozzle 


Ca Cischarge coefficient 


Similarly the second and first stage fiows are given by: 





Quo Ca As > Da (6.8) 
nr 7$ ana P. (6.9) ' 
The desired thrust is given by - 


Š 
U p o.e[. E, C in 


2 | 
(6.10) 





where it is assumed that the load orifice areas are equal to 
the areas of the thire stage power jet anc the factor 0.8 is to 
account for leakage flow out the opposite leg. _ 

From the two stage amplifier, numerical values are 


Known for equations (6.8) and (6.9). “urthermore, since Y, 


is known and assuming a value for cg, equation (6.7) reduces 
to 


253 = (Constant) A 33 (6.11) 


if the numerical values of equations (6.8) ana (6.9) 
are substituted into equation (6.10) along with equation(6.11), 
a quacratic equation L A "S is obtained.  3olving this 
quadrstic equation gives the third stago nozzie area. Since 


all dimensions of the third stage are determined from the 


80 : 





K le width, assuming an aspect ratio for the third stage 


will allow the nozzle width to be determined and hence spec- 


ify the other dimensions oft'the third stage. 


Following the above procedure resulteé in the following 


dimensions for the final design: 





First Stage 
Fower nozzle wiäth = W, = 0.020 in. 
Control nozzle width = 1.Sw, 
Wall setback =  5$.0w| 
Splitter distance = 18, 
Receiver width = Aw] 
Aspect ratio = 1.5w4 


Nozzle erea = 0.0006 ime 


Second Stage 
Power nozzle width = wo = 0.0625 in. 
Control nozzle width = Wo 
Wall setback = 2.5Wo 
Splitter distance = l4Wo 
Receiver width = SWa 
Aspect ratio = 1.5%0wWo 


Nozzle area = 0.00039 in.^ 


Thirà Stage 
Power nozzle width = wç = Ola iaa 
Control nozzle width = Wa 


Wall setback >» 213 
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Splitter distance OW 


Receiver width l.4W 
Aspect ratio 2.084 
Nozzle area 0.0325 in.“ 


( 
A photograph of this device is giown in Figure D.4. 


Ine circular chambers are vortex impedance matchers+ 


modeled efter those used by Compagnuolot®. 

in Figure 6.4 is shown the output thrust of the ampli- 
fler cascade versus supply pressure anc input flow rate. The 
Sucput flow of this final design oscillates aque to the load- 
mE conditions unless the cascace is being driven by an ex- 
ternal signal. 

Eus 6.5 shows the input impecanoe characteristics 
for the right control port when the first stage power stream 


is attached to the right wall. The supply pressure is 1.96 


msi. 


6.5 rersction Jet Cascade Dynamics end Transfer Function 


The reaction jet cascace is characterized by two 
phenomena which ere important when the device is placed in the 
ponerol loop of tne final system: 

1) A hysteresis loop is'present in the switching of 
the thrust from one side to the other 

2) There is a finite time delay between the tine a 
switching signal is received at the control ports until the 
thrust developed is switoheo fron one side to the other. 


The hysteresis loop is due to the characteristics of 
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resection Jet Thrust Versus Supply Pressure 
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bistable devices requiring a finite flow end pressure dif- 


Eu od Lo exist across the control ports to effect switchs 


ing. The half wiäth of the hyster&sis band was determined 





fron ectual measurements with the summer-smplifier set at its 
Peeper supply pressure ang driving the reaction jet cascade. 
The half width was found to be 0.0685 psi. The time celay 
involves two major phenomena: 

1) The time it takes the fluid streams to traverse 
physical distances (transport time) end 

2) The time it takes to charge the interconnecting 
volumes to the required switching levels. 

Actual measurements of this time delay were not mede. 
From theoretical calculations in Appendix B the time delay 
constant was found to bè 0.022 seconds. 


The transfer function of the reaction jet cascade is 


e 


thus given by 


1 


g 


- -Tg 
tT Ze with hysteresis oft 0.0685 
AP psi 
(Seen 
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CHAPTER VII 


CObROL SYSTEM SYNLDIABESIO, J3ru2O0RMÁANOE-AND- &NALYSIS 


W.l öÖvstem Synthesis 





H major part of the problem of synthesizing the system 
components has been accomplished by careful consideration of 
mmemrOlS each ol these components was to play in the final 
system. however, in order to finally bring these components 
together into a working system it is necessary to know how 
each component will behave when actively loaded into the 
cemponent it is to crive in the final design. 

Of critical importance is the setting of the power 
Supply pressure to the Summer-amplifier. This pressure must 
be high enough to prevent the reaction jet cascade from oscil- 
lating and, at the same time, allow for maximum amplirication 
of the rate and position control signals. it can reacily be 
seen that had the reaction jet cascade been cesigned with no 
eonen aeration for this supply pressure setting. that the mags 
nitude or the control signal required for switching the jets 
might have been so large as to preciude the operation of thu 
Systen. 

The rate sensor hac been tested with its outoueS l G 
into passive orifices approximating the control orifices of the 
summer-amplifier. It was thus necessary to recalibrate the 


rate sensor with its outputs,loeded into the summer-amplifier. 
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For accurate results the summer-amplifier must be set at its 
proper supply pressure anc, ,in adaition, must have the second 
Bet oí controls activated.’ ihis is necessary as the input 
we dance of the control nozzles is slightly affectec by the 
EN DIy pressure anc bias ievel at the other controls. ‘The 


steady stste results o? these calibration tests on the rate 


U 


sensor are shown in Figure 7.1. For these tests, the rate 
sensor Output was loaded into the sunner-anplifier with the 
Memwer Supply pressure set at its operating level. The seccnd 
set of controls were not activated. “he results shown in this 
figure were not obtained from an accurate rate table aná are 
considered reasonably accurate only in the range > 0.5 rediens 
per second. 

ln order to set the maximum pressure signal from the 
position sensor”, it was necessary to know: 


) 
1) The gain ratio of the two summer-amplifier control 


ports aná 
2) The final system T/J ratio. 


aS Was pointec out in thes Giscussion of the signal comparator, 
the gain ratio depends on the summer-emplifier power supply 
pressure setting an control bias levels. Thus; to obtain tae 
correct gain retio, a new set of curves Similar to those 


shown in Figure 4.1 were constructed for the proper supply 


pressure setting end control bias levels. 





* The maximum signal occurs at +90 degrees of engular 
rotation. 





EE (Gne ona (20? 





i -1 co (3PY^) 


-6 


Fuure 7.1 
teady State Anrülar Rate Sensor Output 
‘hile Loaded into the Sumner-.rplifier ^ 


! 
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e 


AS Shown in Appendix C.l} the T/J ratio defines the 
switching curve ecuation. “The criteria used to set the maximum 
pressure signal from the position sensor was as follows: If 


! * 
the systom starts with an initial error of os 90 Gegrees and 


e 


O = O radians per second, the maximum rate signal at the 
point of theoretical optimum switching will occur at 6 - 
45 degreeS, je orract switching of the jets, the position 
signal should then equal the rate signal et this point. The 
T/J ratio defines the value of this peak 9 . “rom figure 7.1, 
the pressure signal corresponding to this'value of 8 is deter- 
mined. This then is the value that the position signal must 
 heve at 45 degrees. Knowing the shape of the real switching 
curve and the value of the gain ratio of the summer-amplifier, 
the value of the maximum position signal can readily be deter- 
mined. Í 

In the final system, a pair of adjustable orifices 
were placed between the output of the position sensor and 
the input to the summer-amplifier to make fine adjustments of 
the switching curve. 

Figure 7.2 shows a block diagram of the final control 
system composed of the various transfer functions from the 
previous chapters. A summary of the numerical values for 


the parameters is given below: 


ay 2 1.19 si i 
rad/sec. 


*K nd mp 


- rad 


y 


"These gains have included the gains of the summer- 
amplifier. The summer anplifier multiplied the rete signal 
by six aná the position signal by 4.9 
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To = ome SOC. 
Ta = MOTO CA OC; 
T = 7.1. 10-028 


hysteresis» 0.0688 psi 


J = 9,60 in-oL2-sec^ 


= 


Figure 7.3 shows the! same system in circuit diagram form". 
A photozraph of the complete control system is shown in 


bore 7.4. 
7.2 Svstem Performance 


To obtain performance cata on the system, a rotary 
potentiometer was connected to the control system and the 
output signal displayed on an X-Y recorder. The results are 
shown in Figures 7.5 thru 7.9. 

Figure 7.9 shows tho response of E to a nun- 
ber of position errors. ine reaction jets switched in a near- 
by pernmect manner when responding to counter cloclmise er 
The response to clockwise errors, however, resulted in a 
slight amount of "after end point chattering". The magni- 
tude of the limit is seen to be errsetic, averaging roushly 
G.5 degrees at approximately 0.5 cps. This erratic response 
is due to noise from the proportional amplifiers in the 


system. Considering the T/J ratio of this system, the 


* The symbols used in this figure follow those recom- 
mended by BoothelS. 
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System Circuit Diagram 











Figure 7.4 


Complete Control System and Demonstration Model 
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response time for the errors shown was found to be nearly 
ptimalg 

Pigure 7.6 snows the system operation with a sth, 
NONE Summer-amplitier su ply pressure setting. Sy lowerin 
DE supply pressure, the megnitude oi the hyste:e31s oinG WES 
P oc, resulting in the markeao decresse in the megni tude 
ee limit cycle. Unfortunately, the noise present in Ge 
Eom Coula then switch the reaction jets more readily Te- 
seen in a sligntly less optimel time response to large 
errors because of random switching. 

Sigure 7.7 shows the system performance with the rate 
sensor d'isconnectecd from tre control loop. As pointed out 
more clearly in apperndix C.2, the control system should, 
pusosewudcaliy, Dec completely unstable. because of tren s 
camas in the system, in practite e new limit cycle Was 
reachec at a much higher amplitide and Lower frequency. 

figures 7.8 anc 7.9 snow an interesting mode of oper- 
etione possible with this system. as pointed out in the pres 
cecing chapter, the reaction jet cascace is unstable unless 
Griven, Cue to the heavy loading on che outu c lecs Ol tie 
lust stage. This back pressure is so high that it causes the 
first and second Stage anplifiers also to oscillaee. E 
aue.ne Small control Slows into one side or the first stage 
cause the pulsing to occur more often on one sice than on the 


Other in Clrect proporvics to theme nL uce o hp I o R D 


signal. Thus the system effectively operates in a pulse 
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width moculation mode*. 
ah Fisures 7.9 and 7.9 this pulse wicth moculetion mode 
Wes obtained by lowering the summer-amplifier supply pressure, 
The frequency of oscillation is too high (approximately 
100 eps) for the system to respond to the individusl torgue 
pulses. 
n 


Figure 7.8 shows The response to various disturbances 


while operating in this pulsË wiéth moduletion mode. Of 
@emeorcular note is the disappearance of random ee 
Cue to noise. A certain amount of stiction is noted caused 
by the potentiometer attached to tne system: 

figure 7.9 shows the system response with the rate 

è 

sensor cisconnectea from the control loop. if the load were 
Cc@mp lete ly T OF camping, tne response, cl course, would 
Aca apo Out. a compsrison of this figure with rigure 7.7 
shows ihe functioning of the vortex rate sensor enc the part 
it plays in the stabilization of the sfstem. 
7.5 System snelvsis 

‘This control system falls into the category of non- 
linear control systems because Of the caeng-bang method of 
controlling the load. ‘the most convenient wey of analyzing 
this system, in order to predict the magnituce anc frequenc 
of the limit cycle anc the stability Of the System. 2S 07 


usınp the Cescribi'as-=funetioninemnot. 


Cibson! points out that this approximation methoä is 


"For a more complete explanation of this mode of oper- 
SGlOne See Pere r o pp 


' 
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EN cNecurste in sySvems containing sufficient low pass filte: 
mec, We nature oF en inertial ioac amply fulfills this re- 
quirement es it refuses to respond (filters) to high frequency 
Signals. 

2015S mMetnoo essentially consists Of aSsurimeytheu eMe 
EN Uus to the non-linear element is Sinusoidal. The resulting 
uut is then written in the form of u Fourier series. ‘Then 
it is assumed thst only tne fundamental component of this 
series is important in reletion to the reszonse. The coeficient 
SMEUHIS funcamental frecuency, which is dependent on the mags 
nituce of input signal but indepencent of the frequency, is 


used es an ecuivalent gain (K__) of the non-linear element. 


ec 
The rest of the system being linear is frecuency dependent 
but not dependent on the magnitude or the input signal. By 
clotting the equivalent gein functicn of the non-linear 
element enc the frequency response or Une Trrear elene oo 
ie: plot, a point or intersection may ce found which oce 
Presents a possible point of operation. If the point is 
found to De a staole point, then the frequency and magnitude 
f osocillation cen readily be determined. 

unis metho is applied to the analysis of this system 

in Appendix C.2 the results closely predict both the magni- 


tude aná frequency of the limit cycle shown in Figure 7.5%. 





"ine data on which the calculeticns were based was 
obtained from the system just after cigure 7.5 was recorced. 
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A.l Rate Sensor Design Calculations 





Given Ho = Ow, 


O Š 
Quay = 70 /Sec = 1.221 rad/sec 


AP @ = 0.1 1bf/in* 


max 
SIR an e e 
O = ei lbf -Sec? 
9 ft 


Di 


14 


1.622:107* ft*/sec 


. from equation (2.19) 


| = Mi i a : 


Us > 171.2 ft/sec 






(144) (0.1) 
4(0.00243) (0. 03218) 


From equation (2.21) 


| 2 
; Q = 2 am = tor, | Us 


2 


u 


2 TAS - 4(0.0212) (0.1502 Pi 


g 144 


0.15Y  ft?/sec 9.05 ft?/min. 


e 
H 


103 : 











. 3 h/2 


2 
V Ro 





0150) 0.0581] - 310 
(sera: 7) 00885) 


From reference 3, page 527, we rind that this value of Re is 
off the scale but E is leveling off to roughly 0.66. A 
"factor or safety" was introduced here and E was taken to 


psv... 


Rearranging equation (2.20) gives 


d Ro W $ 

po = 

A US 

~ 2 00.05) 5 5 O 

BUE 

i (090895) (171.2) 

ee e =2 . 

T O = 4.27 dd 10 EU. Spl2 < 10, inc. 


Rearranging equation (2.23) gives 






(Sume Mor P 5 
2 


a = 0604/79 10 
LESS" Cnlculetion of First state 20D) ve ecru ce 


from Figure 2.8, Q = 0.107 in®/sec ané take Q pin 0° 


Co max 








n D P OM ees 
Then V° = GGL 285 in/sec . 


104 








Proz: 


Therefore, VE 2.43 x 10° in^/sec? 





From eouation (2.238) 





a m Ja 
ApV = Qp ECap “p PL 
| 2 
or P Yo P/2 
. p 5 
Cap 


From Corning design data sheets, Cap = 0.70. 


bi = 0.28 psi = 7.7 inches of water. 
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BRREULAT IL UF cho RACION JEF TE LLAT CONSENT 


? 


The time required for the power jet to traverse the cis- 


tance between the nozzle anda the receiver ports is given by 





È = -- i 
' avg 
where 1 <= splitter distance 
Vavg = tne average velocity of a particle 


while traveling between the power 
nozzle exit and the receiver. 


Reference 19 gives the dynamic pressure profile of a sub- 
4 
m C jet et various distsnces from the power nozzle. if 
it is assumed that tne velocity decreases linearly from noz- 


zle exit to the receiver, then Vavg ls given by 


= V. +V a 
Vavg 7 Pisi | 
2 
2 a 
where Vi E D Pg = nozzle exit velocity 
Vs - jet velocity at the entrance to 


the splitter 


The value of UR can be cetermined from the grepns given in 
reference 19 if the splitter distance is known. For the first 
stage: 


Splitter distance = l18w = 0.36 in = 1 
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— TRU A - S 
5 p É Na x 1077) 79.9210 


in/sec 


= & 
P/2 v; | 
š Vig =- v. É Messo in/sco 
Vavg = Vis = Addo 4 108 in/sec 


2 


Lg a 


gs eno” = 73 æ 1076 secs, 





Similarly for the second and third stages: 


Ur 197 £ 107° . secs. 


Ls 
o 


SAT. er SECS. 

Tne time regjuired for the flow*in the control nozzle of 
an intermeciate stage to build to that flow level necessary 
for switching can be approximated from the time constant as- 


Sociated with the volume connecting the two stages ana tne 


resistance of the control port. The time constant is given 
ES 
pe RC 


impedance of the control port lof-sec/in? 


where R 


capacitance due to the interconnecting 
volumes in?/1lof 


O 
i 


The value of R was determined by linearizing the input impeaance 
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curves of each nozzle about an assume operating point. The 


valueWor C is given, for adiabatic conditions, by 


| LO 55e 
PXRT 








where = connecting volume (in®) 


mass density (1bf-sec®/in“®) 


= DO < 
n 


E adiabatic constant = 1.4 


R = gas constant = 2.48 X 10° (in^/sec-94) 


+) 


temperature (?x) ` 


For lines connecting the first and second stage: 


14 


R 


08287 lof-sec/in? 


The value of v was found to equal 0.279 in". sherefore 
! 





c = V 4 08279) 
PERT (1.129:107)(1.4)(2.48X10 
C = 0.013 in” /1of 
Hence 
TEZ .00372 secs. 


The rélationship between the flow into the connecting volume 
(assumed constant) and the flow through the control nozzle 


T 
(variable) is given by 


a RS 
win T S41 


The value Of «in is known fron the nozzle geometry of the 


first stage. «9 is zero before switching occurs. After 
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ening {9 increases until it reaches the value required 


A TUN 


tor switching. The value required for switching was taxen 





to be nearly equal to vin» ‘Therefore, the time required for 
tre flow to Duilà to the switching level is approximately 
NON time constants, or t4 0.0149 seconds. lin a similar 
manner the buildup time between the second and third stage 
was found Tm 0.006 seconds. The rinal time Gelay is 
due to the time recuired for the output flow to traverse 

the length o? the output leg. 3y an elementary calculation 
this was found to equal 636 £ 107° seconds. " 


The total time delav is egual to the sum of all the 


time delays discussed above, 0.022 seconds. 


109 








AMD LL C 


SISIB NeLYo15 


Col Determirine the Zoausticn of the Switchinz Curve 





Consicer a ta ag-bang control system as shown in Figure 
C.l. The egustion of motion describing the operation of this 


system is given by 


9 : sien g] -+ (621) 


ne O is rewritten as 


dt ag at ag 


and substituted for 5 in equation C.l one obtains 


cO aw _ | stsng] (C.2) 
oo i 


Equation 0.2 can now be integrated resulting in 


O: 0,0 ed... Aw 


wi m Torc for  (-) O (C Sa) 

S J 

WF. 29% for ( Q o 
2 d 


= 


These curves represent two families of Dre in the phase 
plane (Wvs e ) where Cy and Ca determine the values of the 
initial conditions. If these curves ere plotted in the phase 
plane, it will be found that the curves corresponding to 


Cj > Co = O pass through the origin. They represent trajectories 


110 








"igure C.1 
Idealized Contactor Servomechanisn Controlling 


an Inertia Load 





Wirure C.2 


Generalized Autonomous Control System Containing 


Y 


a Honlinear Element 


IA 











mh that, if the system is switched at the time one or the 





other curve is intense nm» ontimel reowction of err 3 
meeroduced into tne system will result. 


trus the switching line equations ¿re given by 








wi Ju š I 
Ll. — 2 = Seconc quecrant 
2 switching curve 
e 
NS e 9 = fourth quaärent 


2 Ji Switching curve 


Seo) Limit Cvele analysis 


Ihe block diagram shown in Figure 7.2 can be reducec 
Eure form shown in #igure C.2 ii che following conditions 
exist: 
1j ‘the reference signal is zero 
| 
2) The megnitude of the position signal stays 
within the linear regicn of-the switching curvb." 
Po ue fiture, M represents the nonlineer element in the 
system and G is the combination or eil the linear elements. 


l 


Gibson” shows that amc treen concitien or tue 


existence of sustainec oscillation is 





G( je) z - _1 
Keq 
where Keq = the eouivalent gain of the non- 


lineer element, which is a func- 
ticn of the magnitude of the input 
signal to the ronlinear elezent. 





* This condition is satisfied as shown by the results 
lnw"eeure 7.5. 


112 








&eq for a contactor with hysteresis is given by Gibson as 


. , 
Keo = £ TE) + 3 DIE) * 





where i = Chemeacni twee of tne input to tae 
ecnteetorsrenresenren Ho? A 
and 
g(E) = u 1e (2 E 
° di _ 
b(E)-_ 4T a 
ma E 
E | ' 
where a = the helf wicth ofiysteresisé;&one 


(psi) 


T = the magnitude of the torque output 
(in-ozf) 


For this systen, 


f» 
pt 


0.0686 psi 


By assuming values for — , varíous values of - 1/Xeq 


may be obtained. ror example 





»ssune È = 3 inches of water OTTOS pse 
& = +i l - B. 2 
a Es 
ge. SUD y (0.0688 | ' 
41( 0.1082] 0.1082 


* These values are normalized  Zourier coeficients. 


The term b(E) is needed to take care of phase distortion of 
the fundamental harmonic. 





70 Oi 





8 =: psi 
5 ara s G ın-028 
TER 4" (0.1082) psi 


Keq = 89.5 L 39,49 
-—. : 0.0112 —DSi _ k 39.49 
eq in-ozf 


The linear transfer function G(s) is given by 


k | 
ik SR S zi | = 29 
a(s)= | Kyse” 19, 72, | e773? (gg) 
ToS+1 | 
Js? 


Assume that the time delay of the rate sensor (Ty) is 
zero in order to see what effect rate feedback has on the 


System. “or this case, the linear transfer function becomes 





x DO T: Nu Ky > -TaS E 
G(3): X> y I u u: (0.5) 


J 





S^ (T5841) 


By substituting actual numericsl values gives 


218.59) a | 470: 0225 


” N 
Gie)= (0.1183) (4.89 4,89 
` e 
: S (0.045 +1) 
(ore) 

To show the effect of rate feedback on the stability of this 
system, disconnect the rate feecback (X,=0). For tris case 
equation (C.6) becomes 


GS: 01188 Sas T I 


3° oo) 





Sy substituting S:jG), the values for tne megnituce aac phase 


o. S( WU) as) vàries from zero to inrinity are odteined. I? 


2 


„his 2requency response date is plotted on a poler' plot alon: 


V 


: : š S a ` M = en 

meine megnitude and phase of -i , dB olot Sailer Correia 
ec 

Wei in risure 0.5 will be obtained. lt can be seen that no 





intersection or the two plots occur end hence the system is 
pusocble. | 

Substituting s+iG) into e3gustion (0.6) and plotting thus 
function on a polar plot (Figure 0.4) shows the stebilizing 
effects of rate feecback. It car be seen thet a point of 
Meversection now occurs anc e stedle limit cycle resuits. 
ine effecù of rete fescoack is to introduce phase lead into 
the system which "pulls" the G(jQ@) curve away from the -180 
degrees phase shift line a sufficient am@tunt to siiow this 
intersection to occur. To investigate what effects the time 
mos if) the rate feecdacx will Have on the syster, Soudan on 
(C.4) mey be written in the following form 
A M Deus Ss 


"a en i. ` 
USUS = S æ l + 2e 
J fio =2 








if the rate sensor time delay term is factcrec out of the 
quadratic term anc combined with the reaction jet tire co) np 


“ne result is 


Z a Tan . 
G(3): As E se 5*+ e l e Mtas 
Y . 





(T7 








Pirope C,3 


5 x 
f G(jw) and - E Without Rate !eedbac 
Polam Flot o | 


r 


t 


Im 


figure C,H 


L I 
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UT R ey. yy un eJ ation (0.29), 


be zocec is the increased trensport 


The otner effect to de noted occurs 





U C quacretic term. if tne quedratic terms given in these 
mie ue tions are piottea over the frequency range 0421400, 
it will be “oun that trey differ in magnitude and phase only 
pe very silent amount around tne break Trequency. IDU m us 
memo elect o? the rete sensor time delay is to increase 
@eemearensoort time delay of the over all system. 

“igure u.5 shows a SIS plot of equations (C.5) 09829 


(C.7) where the effect of the rete sensor tire delay cai be, 


secenwson"e clearly. 





Figure C.6 shows polar plots of G(;W) end - E 
The effect of the rate sensor time delay is to ee the 
mMeaenivude of tre limit cycle and decrease the frequency. 
meee this figure, the frequency oF the Limit cyclo Tor tae 


System is precictec to be 9.5 radisas per secoac or l.92 g 


sion for . This procecure gives 


E Te (Q). ) +1 


Z —— 
Qux" E DER 2 x 
ef 222 1(0)% ZI (jW)+1 
Ko Ko S 


, 


where av tre point oz intersection 


: =. = o0. 002 c] 


207 : 





attenuotion eauabions 


(sa2132p) 3Jtut asruj 





7) 


eo 
^e. 


(2.5) ond ( 


Ó 
o 
u) 
cS 
E 
hos 
>> Ze. IAE: m dee | 
We |. .. 4 
NENNEN S | | 1 
o- i o ; | . u 
, ° 
Gà I e e 
m 1 I 


|estt0/0 | Olgor oT apnattiuy 


fnocucney (rıd/sec) 


UN 


~ 
. 


9 


. 
TADO | 


Sinchion G(iw) 


T! m 
LEISTEN 


mer 


v 


“ode Plots of 





Ir 











Magnitude (psi/in-ozf) 
01 222 e DO 





N 


ecuation (©. 7) W) =6 


Ë equation (C,5) 
- 110° “i 
-190° -90° ROS 
120 


B (in, of 1,0) 


| w= 19 
Ly wer ~ 
| Aux |, 
G(jw) CCIW) 


(0) (rad/sec) 


Figure C.6 


Polar Plot of G(jWw) and - — 
m 


L 





ijs, 





lculation gives 


Q = 5.68 degrees 
max 











APPENDIX D 


PHOTOGRAPHS OF COMPONENTS 





Vortex Rate Sensor Showing Details of Sisnal Amplification > 


` 





Figure D.2 Summer-Amplitier 
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Anrular Position Sensor 


t5 ln 





_ 


Picure Ds. 


< 2 


Power Amplifier Cascade and -teaction Jets 





122 








la 


Re 


Se 


4. 


De 


6. 


7, 


9, 


10. 


BIBLIOGRAPHY 


Gibson, J. E., Nonlinear Automatic Control, McGraw-Hill 
Book Company, New York, 1963. 


Taplin, L. B. "Phenomenology of Vortex Flow and its 
Applications to Signal Amplification," Summer Engineering 
Seminars, Pennsylvania State University, Fluid Control 
Systems, July 6-16, 1960. 


Taplin, L. B. and J. F. Hall, "Synthesis of a Pure Fluid 
Flight Control System," AlAA/ION Guidance and Control 
Conference, August 16-18, 1965. 


Sarpkaya, S., "A Theoretical and Experimental Investigation 
of the Vortex-Sink Angular Rate Sensor," Proceedings of 

the Fluid Amplification Symposium, Volume ll, Harry Diamond 
Laboratories, Washington, D.C., October, 1965. 


Dexter, E. M., "An Analog Pure Fluid Amplifier," Symposium 
on Fluid Jet Control Devices, ASME, November 28, 1962. 


Belsterling, Co. å. and K. Tsui, "A Proved Technique for 
Analyzing Proportional Fluid Amplifier Circuits," Çontrol 
Engineering, August, 1965. 


Boothe, W. A., "A Lumped-Parameter Technique for Predicting 
Analog Fluid Amplifier Dynsmics," ISA Transactions, Vol. 4, 
No. 1, January, 1965. 


Katz, S., J. M. Goto and R. J. Dockery, "Experiments in 
Analog Computation With Fluids,” Proceedings of the Fluid 
Amplification Symposium, Volume ll, Harry Diamond 
Laboratories, Washington, D.C., May, 1964. 


Bowles, R. E, and E. M. Dexter, "A Second Generation of 
Fluid System Applications," Proceedings of the Fluid 
Amplification Symposium, Volume 11I, Harry Diamond 
Laboratories, ‘Washington, D.C., October, 1965. 


Blackburn, J. F., G Reethof and J. L. Shearer, Fluid Power 
Control, The M.1.T. Press, Massachusetts Institute of 
Technology, Cambridge, Massachusetts, 1960, 


123 





Bowles Engineering Corporation, "Elements of Pure Fluid 
Systems," Loan Copy from Bowles Engineering Corporation, 
1965. 


Warren, R. W., "Bistable Fluid Amplifiers,” Seminar on 
Fluid Amplification, The Catholic University of America, 
Washington, D.C., June 15-19, 1964. 


Campagnulo, C. J., "A Three Stage Digital Amplifier," 
Proceedings of the Fluid Amplificetion Symposium, Volume I, 
Diamond Ordnance Fuze Laboratories, Washington, D.C., 
October, 1962. 


Campagnulo, C. J. and L. M. Síeracki, "A Digital-Proportional 
Fluid Amplifier for a Missile Control System," Proceedings 

of the Fluid Amplification Symposium, Volume lll, Harry 
Diamond Laboratories, Washington, D.C., October, 1965. 


Hayes, W. F. and C. Kwok, "Impedance Matching in Bistable 
and Proportional Fluid Amplifiers Through the Use of a 
Vortex Vent," Proceedings of the Fluid Amplification 
Symposium, Volume I, Harry Diamond Laboratories, washington, 
D.C., October, 1965. 


Boothe, W. A. and J. N, Shinn,” A CE System of 
Schematic Symbols for Fluid Amplifier Circuitry," Proceedings 
of the Fluid p qaos Symposium, Diamond Ordnance Fuze 
Laboratories, Volume I, Washington, D.C., October, 1962. 


ç 





Warren, R. W., "Pulse Duration Modulation," Proceedings of 
the Fluid Amplification Symposium, Diamond Ordnance Fuze 
Laboratories, Volume 1, Washington, D.C., October, 1962. 


18. Reilly, R. J. and F. A. Moynihan, "Notes on a Proportional 
Fluid Amplifier,” Symposium on Fluid Jet Control Devices, 
ASME, November 28, 1962. 


19.  Peperone, S. J., J. Katz and J. M. Goto, "Gain Analysis of 
the Proportional Fluid Amplifier," Prooeedings of the Fluid 
Amplification Symposium, Volume 1, Diamond Ordnanoe Fuze 
Laboratories, Washington, D.C., Ootober, 1962. 


12% 














D 





